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The first detailed comparison of tlie low-momentum interaction V^ow k with G matrices is pre- 
sented. We use overlaps to measure quantitatively the similarity of shell-model matrix elements for 
different cutoffs and oscillator frequencies. Over a wide range, all sets of Viow k matrix elements 
can be approximately obtained from a universal set by a simple scaling. In an oscillator mean-field 
approach, Viov, k reproduces satisfactorily many features of the single-particle and single-hole spectra 
on closed-shell nuclei, in particular through remarkably good splittings between spin-orbit partners 
on top of harmonic oscillator closures. The main deficiencies of pure two-nucleon interactions are 
associated with binding energies and with the failure to ensure magicity for the extruder-intruder 
closures. Here, calculations including three-nucleon interactions are most needed. Viowk makes it 
possible to define directly a meaningful unperturbed monopole Hamiltonian, for which the inclusion 
of three-nucleon forces is tractable. 

PACS numbers: 21.60.Cs, 21.30.+X, 21.10.-k 



Microscopic nuclear structure studies fall in three cat- 
egories. For local interactions, the Green's Function 
Monte Carlo (GFMC) method [H, leads to exact so- 
lutions of the many-body Schrodinger equation by evalu- 
ation of multi-dimensional integrals in coordinate space. 
The No-Core Shell-Model (NCSM) 0,13 rehes on matrix 
diagonalizations in a harmonic oscillator basis of Nhuj 
excitations with respect to a minimal Ohuj space. Conver- 
gence with Nfko is slow for conventional nucleon-nucleon 
(NN) interactions, which are replaced by effective inter- 
actions that are model-space dependent. Both GFMC 
and converged NCSM methods are limited at present to 
mass number A < 12. The standard Shell-Model (SM) [sj 
restricts diagonalizations to Ohco spaces and treats higher 
excitations in perturbation theory. It bypasses saturation 
problems by using a G matrix [6| calculated at approxi- 
mately the experimental nuclear radius {huj w 40^"-'^/'^) 
and uses experimental single-particle energies. Presently, 
exact SM diagonalizations are possible for all semi-magic 
nuclei, and for A < 70 in full Ohco spaces. 

It has been traditionally assumed that NN interactions 
are strongly repulsive at short distances, and therefore 
require resummations to obtain "pseudopotentials" in a 
given model space. For fifty years the standard in nu- 
clear physics has been the Brueckner-Bethe-Goldstone G 
matrix, which is calculated from a NN potential V by 
summing two-particle ladders outside the model space, 

ri \r Gapkl .-.x 

where are unperturbed (usually kinetic) energies; ij 
and kl denote orbits in the model space and a/3 orbits 
outside it, while the starting energy ojs is treated as a free 
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parameter. Interestingly, the G matrix approach yields 
SM interactions that are, up to an overall multiplicative 
factor of all matrix elements, roughly independent of the 
NN potential and the starting energy used 0, 0) Hi ■ The 
shortcoming of the G matrix is the ill-defined relationship 
between the starting energy and the model space, thus 
precluding ab initio calculations. 

An alternative to the G matrix approach starts 
by noting that conventional NN interactions are well- 
constrained by two-nucleon scattering data only for lab- 
oratory energies Ei^h < 350 MeV. As a consequence, 
details of nuclear forces are not resolved for relative mo- 
menta k > 2.0 fm~^. Starting from a NN potential, 
the high-momentum modes can be integrated out in free 
space using the renormalization group. The resulting 
low-momentum interaction, called V\owk, only has mo- 
mentum components below a cutoff A and evolves with 
it so that all low-energy two-body observables (phase 
shifts and deuteron binding energy) are preserved. For 
A < 2.0 fm^^, all NN potentials that fit the scattering 
data and include the same long-distance pion physics col- 
lapse to a universal Viowfc [§]■ Mowfe defines a new NN 
interaction without a strong core, that can be directly 
used in nuclear structure calculations and therefore elim- 
inates all pseudopotential approximations. 

When only NN interactions are used, all microscopic 
approaches have a common problem, related to poor 
binding and shell formation properties. It reflects in a 
deteriorating agreement with experiment as the number 
of particles increases (active particles for the Ohus SM). 
This leads to the conclusion that three-nucleon (3N) in- 
teractions are necessary. In the case of Viow k , it has been 
shown that chiral 3N forces can be adjusted to remove 
the cutoff dependence and give perturbative contribu- 
tions for A < 2.0 fm^^ in light nuclei In addition, 
these low-momentum 3N forces drive saturation in nu- 
clear matter [l2j . 

In the first part of this paper, we compare Mow fc with G 
matrices. By studying the cutoff and oscillator frequency 
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FIG. 1: (Color online) Top: Correlation plots between Viowfe and G matrix elements in a restricted space of 4 major shells. 
The matrix elements V^^Ju are in MeV for hu = 14MeV and we have highlighted the diagonal elements. Viowfe is derived from 
the Argonne v\g, potential and the G matrix is for Idaho A, computed with a rectangular Pauli operator for 4 major shells 
and starting energy lOs ~ —80 MeV [lC|. Bottom: Correlation plots between Viowfe matrix elements for different cutoffs. To 
facilitate the comparison, we have rescaled the j/-axis set of Viowfe matrix elements by OA^/o-Ay according to the approximate 
scaling law Eq. (Q. The cr-ratios are from left to right: 1.003, 1.061, 1.050 and 1.111. We find even fewer scatters for the T = 1 
matrix elements. All Oab given are O^b- 



dependence of Mow fc matrix elements, we demonstrate a 
tiniversal behavior. A similar behavior exists for G over 
a reasonable range of starting energies. The second part 
is devoted to extracting the Mowfe monopole Hamilto- 
nian, which is then used to calculate binding energies, as 
well as single-particle and single-hole spectra on closed- 
shell nuclei (cs ± 1 spectra) . Many features will turn out 
to be in good agreement with data, and the discrepan- 
cies identify what is expected of 3N forces and how they 
are crucial in heavier systems. Our work shows that the 
mean-field produced by Vjow fe is a valuable first approxi- 
mation that will greatly simplify further perturbative or 
coupled cluster treatment, and the inclusion of 3N forces. 

In Fig. [U we compare V^ow k to G matrix elements in 4 
major shells. We find that both T = and T = 1 matrix 
elements are very similar. For a quantitative cornparison, 
we define the overlaps of interactions A and B 0, S] 

<B = E W^stuA : (2) 

rstuT 

where W/^^ — V/^J^ — Srt ^suW and d2 is the dimension- 
ality of the two-particle space, each state being counted 
[F] = (2 J + 1){2T + 1) times. Here, the interaction V is 
referred to its centroid W, defined by X^rsrl^"] ^rsrs = 0- 
We also introduce normalized overlaps 
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TABLE I: Centroids W , widths a a and overlaps Oab between 
Viow fe matrix elements in 4 major shells for different huj. Viow fe 
is derived from the Argonne wis potential for A = 1.9 fm~^. 
The values of huj correspond approximately from left to right 
to the double-magic nuclei at A = 4, 16, 40, 90, 132 and 208. 
G matrices follow the same behavior. 



with aA — oaa and similarly Ojs ^^r matrix elements 
with the same T . Interactions that differ at most by 
a factor aA/(^B have Oab = 1. The overlaps between 
Viowfe and the G matrix are Oab > 0.99. 

Next, we compare V^ow fe matrix elements for different 
cutoffs. Over the studied range A = 1.3 . . . 3.0fm~^, we 
find again very large Oab overlaps, as shown in the bot- 



torn panels of Fig. [T] To facilitate the comparison, we 
have rescaled the y-axis set of Viow k niatrix elements 
by the widths aA^/a\y. Up to this overall factor, we 
find that Viow k matrix elements are approximately cutoff- 
independent. From TablelH we find a similar behavior for 
sets with different hui at fixed cutoff. These observations 
can be combined in an approximate scaling law 
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(tkjluj U , (4) 



where f//s^„ is a set of two-body matrix elements, ap- 
proximately independent of A and hio. The decrease of 
the overlaps as the range of A or hu increases indicates 
that different parts of the interaction may scale differ- 
ently. However, over a fairly wide range, cfamu follows a 
simple scaling law. 
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Empirically, we have found that the majority of matrix 
elements scale with a = 1 and j3 ~ 1/2. Note that there 
can be cases when Oab is very large, but the scaling 
law is not simple. This happens for instance for different 
starting energies in the G matrices computed from the 
Bonn potential for LUg — ~b . . . — 140 MeV at fixed fiio — 
40 MeV [H, Oab > 0.99 and a ranges from 2.22 to 1.85. 

The preceeding observations suggest that an overall 
multiplicative factor approximately captures the evolu- 
tion of the interaction with A or huj. To this list, one 
would like to add N, the number of hw excitations al- 
lowed in a model space. As shell model calculations with 
any microscopic interaction demand very large N to con- 
verge to the exact result, it is imperative to define ef- 
fective interactions for spaces of lower Nhu) (usually Ohuj 
except for the lightest nuclei). However, for T^owfe at 
A w 2 fm~_\ the particle-particle channel becomes per- 
turbative [l3], the scaling laws in Eq. ^ start operating 
and the Ohuj spaces provide a meaningful first approxi- 
mation. Therefore, we shall proceed by using the bare 
Viowfc for some very simple Ohu; calculations. The ques- 
tion of A^-scaling will be studied in future work, while 
keeping in mind the evidence from G matrices that per- 
turbative renormalizations amount to multiplicative fac- 
tors, though different parts of the interaction may scale 
differently 0,i. 

We start by studying binding energies of closed shell 
nuclei and cs ± 1 spectra. Unless otherwise noted, in 
what follows we use T^ow k derived from the Argonne vis 
potential for A = 1.9 fm~^. We have checked that, for 
small cutoffs, the results are practically independent of 
the precision nuclear force used for Viow k ■ 

As a first approximation the cs ± 1 states are single 
determinants that do not involve configuration mixing 
and are described by the monopole Hamiltonian Hm- In 
the oscillator basis, Hm contains a diagonal and a non- 
diagonal part. The latter is needed for a Haxtree-Fock 
calculation and produces further correlations [5] that will 
be neglected here. 
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FIG. 2: (Color online) Results for the (negative) binding en- 
ergies (BE = (-ffm)) obtained from Eq. (|6]) (Coulomb included 
schematically) by filling lowest oscillator orbits. 



In neutron-proton formalism, the diagonal monopole 
Hamiltonian has a kinetic and a potential part, = 



Hi = R"^ 



ym-r^imSy - Sr^SyS^y) , (6) 



where x,y — n or p, and m^^ is the number of particles in 
orbit r for fluid x. The centroids K-^s^ are defined in 
For an introduction to monopole effects see 14 1 . 

In Fig. [2I we show the binding energies of closed shell 
nuclei. For small Hui, the system is dilute and the interac- 
tion behaves as a contact (5) force, leading to a = 3/2 in 
Eq. We then observe a Hnear dependence on Huj, 

a — I, over the range in which the saturation mini- 
mum should occur once correlations and 3N forces are in- 
cluded. For large hui, the interaction effectively becomes 
long-ranged (a decreases) and the kinetic energy takes 
over. The binding energies are cutoff dependent without 
3N forces. For example, in *°Ca at hw = 12 MeV, we 
have BE/ A = 8.24, 5.89 and 4.52 MeV for A = 1.6, 1.9 
(Fig. [2]) and 2.1 fm~^ respectively. When higher excita- 
tions are allowed, the hui dependence becomes weaker (as 
seen for example in [l^), and an exact result will be in- 
dependent of fku. Fig. [5] indicates that, for medium-mass 
and heavy nuclei, 3N interactions should provide a strong 
repulsion at large Huj. 

The case of ^Li suggests a special behavior in the 
light nuclei in that a minimum is achieved at a rea- 
sonable hu! — 14 MeV, although the total energy is 
still positive due to the kinetic-potential competition: 
{K"^) = 88.0 MeV and (V;^) = -81.3 MeV. Therefore, 
configuration mixing and 3N contributions should lead to 
adding 40 MeV to (F^), a plausible expectation 0, [11]. 
Note that a reduction of Viowfe to the Ofko space would 
demand a factor « 1.4 increase in {V^) to achieve the 
correct result: Again a plausible expectation, suggesting 
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FIG. 3: (Color online) Single-particle spectra on harmonic 
oscillator closures. Lines (meant to guide the eye) join points 
belonging to the same major shell. All energies are measured 
from the largest j subshell in each major shell. We have used 
Eq. ([Sj to refer all Viow k matrix elements to A = 1.9 fm"'^ and 
rescaledby {hw/hcjof = {{p + 2f[{po + 2f + 3]/{po + 2f[{p + 
2)^ +3]^)^ (see text), which corresponds to the mass number 
A appropriate to each major shell with principal quantum 
number p-l- 1. po = 2 (4) correspond to the pf (pfh) shell and 
Hloo = 12 (8) MeV. 



that for medium-mass and heavy nuclei, the 3N interac- 
tions would have to provide an even larger repulsion at 
large hco. 

To obtain further insight, we study the cs ± 1 spectra. 
One known problem of an NN-only description is the fail- 
ure to ensure the N, Z — 28, 50, 82, . . . extruder-intruder 
magicity [B, [IB]- Our results will suggest that this is 
due to two basic shortcomings, the "Z • Z" and "j^ax" 
anomalies, that 3N interactions are expected to remedy. 
These are small effects A^^^ compared to the binding 
energies. 

First, we show the single-particle spectra on top of har- 
monic oscillator closures in Fig. [31 We have scaled all ma- 
trix elements to the fku corresponding to each major shell. 
For this we use a = 2 in Eq. ([5]) to comply with the phys- 
ical imperative that splittings should go asymptotically 
as A'^/^, and we use fuv = 35.59A^^^/{r^)MeV, with 
experimental mean-square radius (r^) = 0.943 A^/'^ (1 -|- 
2/ A) fm^ ti3] and ^ « 2{p + 2 f/3, where p is the princi- 
pal oscillator number at the Fermi level. Our results for 
three values of the cutoff (using the empirical /3 — 1/2) 
and two values of huj are shown in Fig. [3l If the scaling 
laws were perfect, the four patterns should collapse into 
one. They closely do for different cutoffs but less so for 
different hui, indicating that the constant a = 2 is too 
crude. 

The splittings between spin-orbit partners in Fig. [3] 
agree well with experiment, e.g., about 6.0 MeV for both 
the d orbits in ^^O and the / orbits in ^^Ca. The sdg 
and pfh spectra on top of oscillator closures are not di- 
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FIG. 4: (Color online) cs ± 1 spectra in the pf region for a 
wide cutoff range. The experimental 0(?g/2 energy in ^^Ni is 
an estimate. 



rectly available but the necessary information can be ex- 
tracted from the known cs±l spectra up to the Pb region. 
Our results correspond nicely to those determined empir- 
ically [iBl, e.g., 4.6 MeV for the h orbits as in Fig. [3] for 
Hu = 8 MeV, A = 1.9 fm~^. However, what works well 
in is a combination of one-body I ■ s and I ■ I terms. 
The latter is attractive beyond the pf shell, thus favoring 
high I orbits. Fig. [3] shows the opposite: low / orbits are 
always depressed with respect to the I ■ s standard, and 
the effect grows bigger in heavier nuclei. This is the I ■ I 
anomaly. 

Splittings for hole states on harmonic oscillator clo- 
sures are known in ^^O and "^^Ca (« 6 MeV in both). For 
^^O we obtain about half the observed value, a mean-field 
result similar to [l^ . (Note that this hole splitting is not 
the same as the particle one in Fig. [3]). For ■^^Ca we are 
close at 5 MeV, as can be seen in the bottom-left panel 
of Fig. H 

This figure compares the cs ± 1 spectra in the pf shell 
with available data [1^ as summarized in [l^, where 
it is shown that a six-parameter monopole Hamiltonian 
can describe equally well cs ± 1 spectra for A > 60 and 
A < 60. The latter are mostly represented by the pf 
region, which therefore exhibits the basic mechanisms of 
shell formation (see also [1, This one-plus- two-body 
phenomenology reveals what can be explained well by 
realistic NN potentials, and what cannot must therefore 
be ascribed to other (presumably 3N) mechanisms. Our 
task is to establish the distinction. 

The results for low-lying pf levels in the upper panel 
of Fig. |4] are approximately cutoff-independent and the 
minor discrepancies are cured by SM calculations, which 
push up the lp3/2 orbit in '^^Ca, do not change the good 
spectrum in *^Ca, and also push up by 1 MeV the Of 5/2 
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orbit in ^^Ni 0, For the high- lying sdg shell levels, 
the splittings between spin-orbit partners agree with the 
empirical values, but the I ■ I anomaly is much stronger 
than in Fig. [31 and partly responsible for the ^ 10 MeV 
underbinding of the 0179/2 particle orbit in ^^Ni. A sim- 
ilar discrepancy shows for the O/7/2 hole orbit in ''^Ca, 
underbound by ^ 10 MeV with respect to its sd part- 
ners. Such shortcomings are responsible for the failure of 
NN-only interactions to ensure the N, Z ~ 28, 50, 82 . . . 
extruder-intruder closures. This can be directly checked 
through the standard measures of magicity, the gaps de- 
fined as g{cs, x) — 2BE{cs) — BE{cs + x) — BE{cs — x) 
(x = n or p). The ground state spins are those of Fig. |4l 
Rounded calculated (experimental) gaps are (in MeV): 
.g(4°Ca,n) = 13(7), 5(^^Ca,p) = 17(6), 5(^^Ca,n) = 
-0.4(5) and gi^'^Ni, n) = 0.5(6). It follows that extruder- 
intruder closures are non-existent, and harmonic oscilla- 
tor closures are too strong. A related problem is that 
the 0^5/2 hole orbits in "'^Ca and '^^K are underbound by 
^ 4 MeV with respect to their sd counterparts. On the 
contrary the lsi/2 — 0^3/2 splittings are quite good: Al- 
though too large by some 1.5 MeV in "^^Ca, in "'^Ca the 
splitting is drastically reduced, close to what the data 
demand, and further reduced in ^^K, now very close to 
experiment. 

Thus, we find that when the largest j orbit in a major 
shell fills, it binds itself and contributes to the binding of 
the largest j orbits in neighboring shells in a way that NN 
forces fail to reproduce. This is the jmax anomaly. The 
necessary intra-shell self-binding to cure it is now well 
understood in terms of a 3N mechanism [2^ . A mecha- 
nism to resolve the cross-shell binding problem detected 



here remains to be found. 

In summary, we have shown that Viow k and G matrix 
elements are quantitatively similar, but Viowfe as a free- 
space potential is far easier to use in many-body calcula- 
tions. As a consequence of the similarity, it is possible to 
build on the successes of the G matrix approach, without 
the drawbacks due to the ill-defined starting energies and 
other limitations on the SM by hard potentials. Viowfc 
leads to matrix elements that are approximately cutoff 
and oscillator frequency independent up to an overall 
scaling with the width of the interaction. The scaling 
properties associated to the universal behavior of Viow fc 
have been tested in cs ± 1 spectra. The soft nature of 
^ow k allows direct monopole estimates of binding ener- 
gies and cs ± 1 spectra, and in an NN-only description, 
Viowk reproduces many features of the cs ± 1 spectra. 
Moreover, our results suggest that, apart from satura- 
tion, the main problem that demands a 3N interaction is 
related to extruder-intruder shell formation. 

Finally, it is worth mentioning that ongoing NCSM 
calculations ^2§\ nicely confirm the validity of our Ohuj 
explorations. In particular, the saturation patterns for 
A = 6, 12, 16 follow those of Fig. [2] within an overall en- 
ergy scafing and the spectra of ^ = 15, 17 are the same 
within few hundred keV. They provide some first sugges- 
tions on how nuclear spectroscopy can help to constrain 
the 3N interactions that prove crucial for a microscopic 
understanding of nuclear many-body systems. 
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